Spin squeezing and Schrodinger cat generation in atomic samples with Rydberg 

blockade 



Tomas Opatrny 1 and Klaus M0lmer 2 
1 Optics Department, Faculty of Science, Palacky University, 
17. Listopadu 12, 77146 Olomouc, Czech Republic 
2 Lundbeck Foundation Theoretical Center for Quantum System Research, 
Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark 

(Dated: March 6, 2013) 

A scheme is proposed to prepare squeezed states and Schrodinger cat-like states of the collective 
spin degrees of freedom associated with a pair of ground states in an atomic ensemble. The scheme 
uses an effective Jaynes-Cummings interaction which can be provided by excitation of the atoms to 
Rydberg states and an effective J x interaction implemented by a resonant Raman coupling between 
the atomic ground states. Both dynamical evolution with a constant Hamiltonian and with adiabatic 
variation of the two interaction terms are studied. We show that by the application of further 
resonant laser fields, we can suppress non-adiabatic transfer under the time varying Hamiltonian 
and significantly speed up the evolution towards a maximally squeezed, J z — 0, collective spin state. 
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I. INTRODUCTION 

Atomic ensembles hold great potential for precision 
metrology [U [3] as well as for quantum information pro- 
cessing [3H5], where one can take advantage of the dis- 
crete character of spin states of an individual atom as well 
as of the quasi-continuous state structure of large atomic 
samples. In metrology, the precession of the atomic spin 
is used for time measurement in atomic frequency stan- 
dards and for high precision probing of, e.g., magnetic 
and electric fields. Engineering collective spin squeezed 
and entangled quantum states makes it possible to sup- 
press noise and increase the measurement precision and 
sensitivity 1,2.6 9J . The ladder of quantum states of the 
collective spin is formally compatible with the structure 
of quantized radiation modes, and the conventional atom- 
light interaction thus enables the construction of matter- 
light interfaces for quantum communication, where one 
may also benefit from ensembles being initially prepared 
in squeezed and entangled states [10] . 

Atomic ensembles are readily manipulated with laser 
fields that drive transitions between internal atomic 
states. To produce states with entanglement between 
different atoms, however, more sophisticated interaction 
schemes are necessary. To prepare such states, one ei- 
ther needs to bring the "nonclassicality" from outside, 
e.g., by absorption of nonclassical states of light [TT] or 
by the quantum back-action of an optical measurement 
of a collective atomic observable [7J [9] , or one may in- 
duce suitable interatomic interactions within the ensem- 
ble. The former approaches have their limits as not all 
quantum states of the light field are easily available and 
the measurement sensitivity has to beat the usual shot 
noise limit; and to have sufficient atom- light coupling 
strength, they typically work only for very large sam- 
ples. The latter approach is not easy because atomic en- 
sembles are typically dilute and the atoms move around 
randomly: the interactions among the atoms, needed for 



quantum state preparation, are thus also a source of de- 
coherence. 

In this paper we propose to use the interaction be- 
tween Rydberg excited states [T^] to correlate the atoms 
within an ensemble in a controllable way. There have 
been various proposals to use the Rydberg interaction 
for quantum state manipulations and for quantum infor- 
mation processing |13H21) . As in a number of previous 
proposals, we make use of the so-called blockade effect, 
where the excitation of one single atom shifts the energy 
and thus prevents the resonant excitation of other nearby 
atoms. This effect on a single pair of atoms has been ob- 
served and used to create entangled states and perform 
quantum gates [22J H3| and in large atomic samples it 
leads to a significant suppression of the excitation num- 
ber and number fluctuations under resonant irradiation 
[25] . We assume an ensemble of atoms confined to 
a spatial volume so that all atoms are within the inter- 
action range of the others. The quantized occupation 
of two stable atomic ground states is then represented 
by harmonic oscillator degrees of freedom while the re- 
striction of the number of Rydberg excited atoms to zero 
and unity permits a mapping of this degree of freedom 
on an effective two-level system. Moreover, the coher- 
ent coupling of an atomic Rydberg and ground state has 
the form of the Jaynes Cummings (JC) model [55] of the 
interaction between a two-level atom and the quantized 
radiation field in quantum optics. 

We will investigate to what extent the well known the- 
oretical potential to generate squeezed states [37J HE] and 
Schrodinger cat states [H] of light applies also to the case 
of collective atomic spins. Apart from the formal anal- 
ogy, which allows us to directly apply the results of the 
JC coupling, our scheme takes advantage of some special 
properties of the atomic oscillator system as compared to 
the quantized field: (i) the atoms populate two ground 
states and they can be prepared in the spin-coherent state 
at the equator of the Bloch sphere where the binomial 
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FIG. 1: Atomic level scheme: two low lying states |o) and \b) 
are coherently coupled by resonant Raman laser fields with an 
effective Rabi frequency Q, and they are separately coupled 
to two Rydberg states 1 7*1,2} with effective Rabi frequencies 
^1,2- The dots represent a possible population of the atomic 
levels. 



number statistics has smaller fluctuation than the corre- 
sponding coherent state of the harmonic oscillator with 
the same mean occupation, (ii) the JC interaction can 
be switched on and off and both the intensity and phase 
of the JC Hamiltonian can be relatively easily varied in 
time, (Hi) one can couple each of the low lying atomic 
levels to separate Rydberg state and thus eliminate some 
unwanted asymmetries generated in the standard JC sce- 
nario, and (iv) one can simultaneously act by resonant 
Raman fields to adjust the atomic ground state. There 
are of course numerous further practical differences be- 
tween squeezing of light and atoms, both concerning the 
lifetime of the states produced and their possible appli- 
cations. 

The paper is organized as follows. In Sec. [TTJ we intro- 
duce the physical system and our notation. In Sec. |III| 
we present numerical and analytical results for dynamical 
squeezing with the effective JC interaction 
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In Sec 

we determine the eigenstates of a Hamiltonian with both 
Raman and JC coupling terms, we show that it adiabat- 
ically connects a spin coherent state with a maximally 
spin-squeezed state, and we show that application of 
judiciously chosen additional couplings suppresses non- 
adiabatic processes and permits rapid evolution along 
the desired eigenstates. In Sec. [V] we analyze the pro- 
duction of Schrodinger cat states by the JC interaction. 
In Sec. IVIl we summarize our results. 



II. SYSTEM DESCRIPTION 



states I a) and \r-y) and between \b) and \r%) are shown in 
Fig. [I] we will also make use of coupling between states 
I a) and |r 2 ) and between \b) and |ri). There are N atoms 
in the sample and we assume that their interaction pa- 
rameters with the fields are identical, so that if initially 
all the atoms are prepared in the same ground state, say 
I a), any action of the fields will lead to states confined 
to the manifold of symmetric superpositions of all the 
individual atomic states. For co-propagating laser fields, 
the propagation phases over the spatial extent of the en- 
semble play no role for the ground state coherence, while 
they may be absorbed in the definition of the the opti- 
cally excited states to render the field couplings and state 
amplitudes fully symmetric under permutation of atoms. 
Let us denote by \n a , n^, n rl , n r2 ) the symmetric state in 
which n a atoms are in state |o), rib are in state \b), and 
n r i, r 2 atoms are in the respective Rydberg state 1^1,2)- 

A single atom of the sample can be laser-excited to 
a specifically chosen Rydberg state with large principal 
quantum number n[T2]. Due to the strong dipolar cou- 
pling between two Rydberg excited atoms, excitation of 
another atom to the same Rydberg state would not be 
resonant with the same laser field, but would require a 
detuning given by the dipole interaction energy. The van 
der Waals interaction gives rise to frequency shifts of few 
kHz for n — 45 Rydberg excited Rb atoms separated by 
few micrometer distance, while Forster resonances occur 
for special principal quantum numbers and leads to shifts 
up to few MHz, e.g., for Rb atoms excited to 40p3/2-states 
up to 10 micrometers apart [21]. This shift is at the root 
of the Rydberg blockade mechanism: when an ensemble 
enclosed within a few micon sized volume is uniformly ex- 
cited with a laser field resonant with a specific Rydberg 
state, only one excitation occurs within the ensemble, 
and as long as the excitation Rabi frequency is well be- 
low the blockade shift, further excitation is energetically 
blocked. Our Hilbert space is thus restricted to states 
\n a , rib, n r i, n r 2) where n r \. r i are confined to and 1 and 
where n a + rib + n r \ + n r2 = N. As can be checked, this 
restricted Hilbert space has dimension 47V. 

The dynamics between states |o) and \b) can be de- 
scribed by means of the collective spin operator J = 
^ fc Sk, where the index k denotes individual atoms and 
the individual atomic spin components are 



We consider an ensemble of atoms with the level struc- 1 

ture shown in Fig. [I] two lower states \a) and |6) are S x = ^ 

(|o)<6| + |6)<o|), (1) 

coherently coupled to two Rydberg states |ri,2) by laser j 

fields (in most experiments coupling to Rydberg states is S y = ^ ( — \ a )fi\ + 1^) \ a l)' (*) 

achieved by two photon transitions via an intermediate j 

excited atomic state), as well as to each other by a co- S z = -(|a)(a| — |6)(6|). (3) 
herent Raman coupling. We assume that the laser fields 
at different frequencies are derived from a single master 
or that their phases are otherwise locked to a common 

reference. Although only the Rydberg coupling between The collective spin components can be represented by 
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means of creation and annihilation operators as 

J x = ^(otft + oftt), (4) 

J y = l -(a?b-atf), (5) 

J z = \{<Ja-b%), (6) 

where the operator a destroys an atom in state \a), a* 
creates an atom in state |a), (and similarly for 6), and 
the commutation relations are [a,<v] — [b,b<] = 1. 

Hamiltonian terms proportional to the components of 
J can be realized by electromagnetic fields near resonance 
with the a — b transition: the detuning of the frequency of 
the field then gives the J z contribution and its amplitude 
£1 and phase ip determine the J XiV contributions to the 
Hamiltonian. 

The single excitation allowed in either of the Rydberg 
states is distributed in a symmetric fashion as a superpo- 
sition state with the same excitation amplitude on each 
individual atom, and the coupling of the ground state 
with n atoms to all the components in the excited super- 
position state leads to an enhancement by the factor yfn 
compared to the single atom coupling. As the Rydberg 
excitation reduces the ground state atom number by one, 
this process couples states of the form \n, 0) and \n— 1, 1), 
and the y/n dependence of the coupling strength is in 
exact agreement with the matrix element of the Jaynes 
Cummings Hamiltonian Hjc = ga + <j- + g*a<j-. The 
Rydberg excitation of atoms subject to the excitation 
blockade thus implements the JC Hamiltonian, where the 
oscillator degree of freedom represents the ground state 
population and the two-level system represents the Ry- 
dberg state population. In our model, the Hamiltonians 
that couple the two lower states and the Rydberg states 
can be described as 

Hjci = fiio^ +^* 1 a t cri 1) , (7) 
Hj C 2 = £l 2 baf +ft*foVi 2) , (8) 

where the operators o± create or destroy an atom in 

(2) 

the Rydberg state \ri), and similarly for a± . If the 
Rydberg coupling field £lj is detuned from resonance by 
5j , the Hamiltonian will also contain a term proportional 

to SjaP where = af >a U} - a^exf. 

III. DYNAMIC SQUEEZING 

As discussed in detail in [35], the JC interaction with 
a two-level system can evolve coherent states of a har- 
monic oscillator into amplitude squeezed states. This is 
most easily understood if we assume that the two-level 
system occupies one of the dressed state superpositions, 
(|0) ± \l))/y/2, which diagonalize the interaction with a 
classical resonant field with a real Rabi coupling. The 
time evolution of a product state of such a dressed state 



and an initial coherent state of the oscillator is approx- 
imately given by a simple phase evolution of each os- 
cillator eigenstate component, corresponding to the Rabi 
interaction energy, which is proportional to i/n. This im- 
plies a rotation and a twisting of the coherent state am- 
plitude distribution in phase space towards a squeezed 
state. Unlike the evolution under the usual squeezing 
Hamiltonian which is quadratic in the creation and an- 
nihilation operators, the squeezing in the JC system is 
not ideal, and the circular phase space distribution does 
not evolve into an ellipse, but rather into a banana-like 
shape [25] , 

To produce a spin-squeezed states in an analogous 
manner, we initialize the atomic sample with all atoms in 
an even superposition of states \a) and |6), corresponding 
to a spin coherent state with maximum eigenvalue of the 
J x operator. Thus, while each atom is in the superposi- 
tion 2 -1 / 2 (|a) + \b)), the collective state is 



hAo)=2- w / 2 ]T J{ j|n a ,A-n a ,0,0). (9) 

n a — V ^ a ' 

Next we want, for all values of n a and rib, within the 
corresponding binomial distribution, to prepare the two- 
level dressed states with states \r\) and \r 2 ) populated 
and unpopulated with equal amplitude. Within the quite 
narrow distribution of n a n,\, tt/2 excitation pulses on the 
two Rydberg transitions, may accomplish this state with 
adequate precision. The precision may be further en- 
hanced if the Rydberg state is excited by an adiabatic 
chirp of the Rydberg exciting laser detuning towards res- 
onance. During the chirp, phase factors depending on 
n a n,) are accumulated, but since their values are known 
and since they contribute in the same manner as the 
phases we shall need for the squeezing process, we shall 
assume that they present no problem for the correct ini- 
tialization of the system. 

After having prepared the initial state, the system 
evolves with the resonant Jaynes-Cummings Hamiltonian 
Hjc = Hjci + Hjc2, where we choose fli — fl 2 = £ljc- 
The system thus corresponds to a pair of standard JC 
systems in which the initial state of each system approx- 
imately corresponds to the "field" in a coherent state 
with the mean number of photons (N — l)/2 and the 
"atom" in the superposition of the ground and excited 
states. Our system differs from the standard JC by 
two aspects: (i) the "photon" numbers n a> h and the 
"atom excitation" numbers n r \^ r 2 are anticorrelated by 
n a + rib + n ri + n r 2 = N, and (ii) having n r i^ r 2 fixed, 
n a ^b have a binomial rather than a Poissonian distribu- 
tion. The anticorrelation is very important for the shape 
of the resulting state: whereas the a-mode component 
with n a = n + An atoms acquires a phase proportional 
to \/n + An ss + ^ — f §=V), the 6-mode component 
with rib = fi — An atoms acquires a phase proportional 
to \Jfi — An w y/fl(l — — ^=V) (assuming the mean 
excitation number ft — (N — l)/2). These two contri- 
butions add so that the terms linear in An, as well as 
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FIG. 2: (Color online) Probability of finding n a atoms (wide 
blue bars) in state \a) in a sample of N = 64 atoms that 
underwent the dynamical squeezing procedure. The narrow 
red bars represent statistics of the spin coherent state, for 
comparison. 
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FIG. 3: (Color online) Q-function of the resulting squeezed 
state in a sample of N = 64 atoms after the dynamical squeez- 
ing procedure. 



all other odd power terms, cancel. Apart from an unim- 
portant constant, one is left with a phase evolution pro- 
portional to the even power terms of An, in particular, 
4>k -Q, Jc tAn 2 /(4n 3 / 2 )+0(An 4 ). If terms beyond third 
order can be neglected (for An <g; n) , the evolution in the 
a and b modes corresponds to that caused by a Hamil- 
tonian oc J z that leads to spin squeezing by "one-axis 
twisting", as shown in [6\. Thus, whereas in the stan- 
dard JC model the coherent state rotates in the phase 
space (due to the linear term in An) and aquires a "ba- 
nana" shape during the squeezing procedure, in our sys- 
tem states starting at the equator of the Poincare sphere 
do not rotate and they stay symmetric. Moreover, in 
comparison with the coherent state with the standard 
deviation of photon number = \ffi, the binomial distri- 
bution is narrower with the standard deviation = yJh/2. 
This also contributes to render the higher powers of An 
in the evolution less important than in the standard JC 
model, and we reach better squeezing than in the stan- 
dard JC model with the same initial mean excitation 
number. 

The time for reaching the maximum squeezing can be 
estimated by considering that states with the excitation 
number around n + An rotate in the phase space by a 
phase angle ~ |f2jc|iAn/(2n 3 / 2 ) and the initially spin co- 
herent state becomes maximally squeezed when the con- 
tributions from the ends of the distribution n ± y/n/2 
have rotated to the opposite angles ±7r/2. Therefore the 
time of maximum squeezing is ~ iV/|f2jc|- This esti- 
mate is well confirmed by our numerical simulations of 
the process. 

The combined JC operators induce a nonzero correla- 
tion (JyJ z + JzJy) equivalent to a tilted uncertainty el- 
lipse which can be rotated by means of a Raman transfer 
Hamiltonian cx J x , so that the two ground state popula- 
tions acquire subbinomial distributions. The state of the 
ensemble contains contributions of the Rydberg states, 



which can be transferred to the ground states by fre- 
quency chirping the Rydberg coupling fields out of reso- 
nance. The state is now squeezed in the J z variable, i.e., 
it has suppressed fluctuations in the difference of atom 
numbers n Q — rib i n levels a and b. 

In Fig. [2] we show results of our numerical simulation of 
the procedure for squeezing a state with N = 64 atoms. 
The final fluctuation is AJ Z — 1.08 which corresponds to 
11 dB squeezing of noise below the fluctuations of a spin 
coherent state having AJ Z = 4. In Fig. [3] the Q-function 
of the resulting state is displayed (the Q-function is the 
trace of the product of the density matrix of the state 
with the density matrix of the spin coherent state). 



IV. ADIABATIC SQUEEZING 
A. Hamiltonian spectrum 

For an even atom number, TV, it is in principle possi- 
ble to reach a collective J z = eigenstate with A J z = 0, 
i.e., n a = rib, whereas for odd N the maximally squeezed 
state has A J z = 1/2 associated with the two equally pop- 
ulated components n a = n&±l. To assess such maximally 
squeezed states we consider adiabatic transition from an 
extremal eigenstate of operator J x (i.e., a spin coherent 
state) to an extremal eigenstate of Hjc- Since both of 
these operators and any weighted combination of them 
can be implemented with suitable laser fields, we are in 
principle able to drive this adiabatic transition, and let 
us show that the extremal eigenstates of Hjc are indeed 
maximally squeezed. 

The eigenstates of Hjc are dressed states of a form 
analogous to the standard JC model. Assuming a real 
coupling strength Vljc, the usual pair of eigenstates gen- 
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eralizes to four states, 



|<+ nb) ) 



|V>(_"°j nb) ) 



\^X nb) ) 



1 



(|n„,n 6 ,0,0) + |n a -l,n 6 ,l,0) (10) 



|n a ,n 6 - 1,0, 1) + \n a - l,n b - 1, 1, 1)) , 
i(|n a ,n 6 ,0,0) + |n o -l,n 6 ,l,0) (11) 
|n a ,n b - 1,0, 1) - \n a -l,n b - 1, 1, 1)) , 
i(|n a ,n h ,0,0)-K-l,nb,l,0) (12) 
|n a ,n 6 - 1,0, 1) - \n a -l,n b - 1, 1, 1)) , 
i (|n o ,n&,0,0) - |n a - l,n&, 1,0) (13) 



- |n a , rib - 1,0, 1) + \n a — l,n b — 1, 1, 1)) 
with the energies 



(n a ,n b ) 
+,+ 

(™a,«b) 
+ ,- 

(™a,™b) 
-,+ 

(™a,™b) 



(14) 
(15) 
(16) 
(17) 



For even N the extremal eigenstates have ener- 
gies ±QjcV2N and correspond to n a — n b — 
N/2, while for odd N each of the extremal states 
is doubly degenerate, e.g., the states with max- 
imum energy Q rc (y/(N + l)/2 + y/(N - l)/2) are 
| ^«iv + i)/2,(iV-i)/2) ) and ^((iv-ijA^+D/a)^ Thug; in 

the extremal energy eigenstates, the populations of lev- 
els a and b are as close to each other as possible. Eigen- 
values of the combined Hamiltonian xHjq + (1 — x) J x 
for x between and 1 and N = 16 are shown in Fig. [4] 
As can be seen, an extreme eigenvalue of J x is smoothly 
transformed into the corresponding extreme eigenvalue 
of Hjc when the parameter x changes between zero and 
unity. Thus, by properly choosing functions f 1,2ft) one 
can steer the state by adiabatically changing the Hamil- 
tonian fi(t)J x + f2(t)Hjc such that an initial spin coher- 
ent state (i.e., extremal eigenstate of J x ) evolves into an 
Hjc extremal eigenstate. After the adiabatic transfer, 
one can get rid of the Rydberg excitations by sweeping 
the Rydberg coupling helds out of resonance, and if N is 
even the final state is \N/2, N/2, 0, 0) . 

If N is odd, the final superposition state, (\(N + 
l)/2, (N - l)/2, 0, 0) + \(N - l)/2, (N + l)/2, 0, 0))/V2, 
carries a phase that depends on the orientation of the 
initial phase coherent state. With imballanced JC 
strengths, one may in principle prepare one of the states, 
e.g., J z = 1/2. 



B. Compensation for nonadiabatic transitions 

The adiabatic scenario holds the promise to generate 
much better squeezing than the dynamic time evolution 
with the fixed JC interaction, but in the vicinity of the 
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FIG. 4: (Color online) Eigenvalues of Hamiltonian xHjc + 
(1 - x)J x for TV = 16. 



extremal states, Hjc has a very narrow level spacing 
(oc QjcN~ 3 / 2 ). Hence, too fast parameter changes will 
couple different adiabatic eigenstates and thus decrease 
the resulting squeezing. One can deal with this prob- 
lem by a careful choice of the functions ftpty) such that 
the trade-off between transitions to unwanted states and 
the speed of the process is optimized. To maintain adia- 
baticity even for moderate atom numbers will, however, 
be incompatible with the finite Rydberg state lifetime. 
We propose here a strategy to actively compensate for 
the nonadiabatic transitions by additional time-varying 
Hamiltonian terms. For a given time varying Hamilto- 
nian H (t), the well known non-adiabatic coupling terms 
due to the time dependence of the adiabatic eigenstates 
are of a form fully equivalent to the application of an ex- 
tra Hamiltonian in a fixed basis, and as proposed in [30j . 
they can therefore be canceled by applying precisely the 
negative of that Hamiltonian, which is explicitly given as 



Hi(t) 



\m)(m\d t H \n)(n\ 



E n — E 



(18) 



A system subject to the combined Hamiltonian Ho(t) + 
Hi (t) evolves exactly along the instantaneous eigenstates 
of H a (t) . It is of course in general not easy to provide ex- 
actly the Hamiltonian H 1 (t) from physically available in- 
teractions. Even though one may not be able to construct 
the full Hamiltonian H\(t), one can, however, attempt to 
preserve as close as possible the evolution of the extremal 
adiabatic eigenstate \i/)o(i)) of Ho(t) by the application of 
a judiciously chosen pertubation of the system. Among 
the operators that can be naturally implemented, and 
which couple the adiabatic states so that they can be 
used to counter the non-adiabatic transitions, are 



H 



H (y) 
11 JC 

(y, cross) 
JC 



(a*W a i a W + ba f _ )(19) 
[ba^-b^^+aa^-a^^). (20) 
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Note that Hjq has a similar form as Hjc (uses the same 
transitions between the lower levels and the Rydberg 
states), but the phases of the coupling fields are differ- 
ent. The operator jj^ lass ^ uses fields driving "cross" - 
transitions between states a and r 2 and between states 
b and r 2 (not shown in Fig. [I]). To partially compensate 
for the nonadiabatic transitions, we thus suggest to drive 
the system with the Hamiltonian 



H(t)=H {t) 



ai {t)H% 



a 2{t)H JC 



where 



H (t) = h(t)J x + f 2 (t)Hj C 



(21) 



(22) 



is the original Hamiltonian and the functions oti l2 are 

chosen such that the norm of the vector [ai{t)H^ + 

a 2 (t)Hj V c class ^ — Hi(t))\ipo(t)} is minimized. This results 
in an explicit formula for ax,2 involving only matrix ele- 

iv) Tr(v,cross) 



ments of the operators H y J( j, H JC 



and Hi{t) in the 

state \ipo{t)) (see Appendix |A|) 

We note that the availability of a number of different 
interaction Hamiltonians with variable strengths is the 
prerequisite of optimal control theory [3TJ [32] , and that 
strong numerical methods exist to identify the fastest 
and most reliable route towards desired final states. Our 
approach towards a useful choice of parameters is almost 
with certainty not the optimal one. We believe, however, 
that it offers an interesting, fast and explicit protocol, 
and that it retains a physical interpretation, which guides 
our efforts to choose the few, most relevant interaction 
terms. 



C. Results 

The procedure and its results are shown in Figs. [5J- 
for a system with N = 15 atoms. The rate of the tran- 
sition increases with the magnitude of parameters /i,2- 
The fastest process would occur if one switches the JC 
coupling f 2 (t)Hjc to the maximum possible value and 
then slowly turns off the J x part of the Hamiltonian, the 
rate of turning J x off increases with the magnitude of f 2 . 
Since the JC coupling to the Rydberg state cannot be ar- 
bitrarily strong (because one may, e.g., excite more than 
one atom to the Rydberg state if fijc becomes compara- 
ble to the blockade splitting), it is natural to choose the 
maximum allowed value of the JC coupling as the princi- 
pal restriction in the optimization procedure. In our very 
simple scenario, the algorithm in each time step chooses 
the change of parameters fi j2 such as to approach the tar- 
get values (fi — >■ 0, f 2 1, see Fig. [5]) as fast as possible, 
while having the amplitudes of nonadiabatic transitions 
within a predetermined tolerance interval. 

The resulting time dependent spectrum of the Hamil- 
tonian Ho(t) is shown in Fig.[6j The slow approach to the 
terminal state is necessary due to the tight level spacing 
of Hjc- The figure also shows as a green (red) curve the 
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FIG. 5: (Color online) Evolution of parameters fi(t) and /2(f) 
in the Hamiltonian Ho(t) chosen such as to minimize nona- 
diabatic transitions provided the maximum JC coupling is 
limited and the state should proceed to the final stage as fast 
as possible. The time is given in units l/|f2jc| correspond- 
ing to the maximum JC coupling. The number of atoms is 
TV = 15. 
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FIG. 6: (Color online) Time evolution of the spectrum of 
Hamiltonian Ho (t) with parameters as in Fig[5] The red curve 
(close to top, marked with "x" in detail) is the mean value of 
energy calculated during the evolution. For comparison, the 
green curve (marked with dots in detail) is the mean value of 
energy for the procedure without the compensation of nona- 
diabatic transitions. 



mean value (Ho(t)) in the evolving state without (with) 
compensation for the nonadiabatic transitions according 
to Eq. (21 ) . The perfect agreement of the red curve and 



the extremal eigenvalue of Ho(t) demonstrates that our 
compensation significantly improves the adiabatic follow- 
ing. 

In Fig. [7]we show the time dependence of the compen- 
sation parameters ai t2 of Eq. (21 1. It turns out that the 
available operators only have substantial overlap with the 
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FIG. 7: (Color online) Time dependence of parameters 0:1,2 
of the Hamiltonian Eq. ( |21[ ) chosen such as to minimize the 
nonadiabatic transitions. 
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FIG. 8: (Color online) Time evolution of the squeezing param- 
eter S in the dynamic regime and in the adiabatic procedure 
with compensation of non-adiabatic transitions, cf., Eq. (21 1. 



Berry compensation Hamiltonian H x of Eq. pi in the switched on for time r, where 



initial stage when Hq is dominated by J x . In the later 
stages this overlap becomes very small and the compen- 
sation becomes virtually inefficient. However, even with 
such a limited option for compensation the additional 
Hamiltonians enable us to perform the transition sub- 
stantially faster, by a factor of ~5-10. 

We may terminate the process described in Figs. |4j [7] 
at any instant and apply a chirp of the JC interactions to 
remove the Rydberg excitation. The resulting time de- 
pendent value of the squeezing parameter S = (A Jf ) /N 
is shown in Fig. [8] Note that a phase coherent state with 
(J z ) = (i.e., our initial state) has S = 1/4 correspond- 
ing to the binomial distribution, and any state having 
S < 1/4 is squeezed. As can be seen, the dynamical 
procedure squeezes the state faster but after reaching a 
certain minimum value the parameter S returns to val- 
ues of very poor squeezing. The adiabatic process works 
more slowly but leads to a much deeper squeezing. 



V. SCHRODINGER CAT GENERATION 

To generate superpositions of coherent spin states with 
opposite spin directions, one starts with the same initial 
state as in the preceding sections, but uses coupling only 
to one of the Rydberg levels. The scheme has some simi- 
larities to that of standard JC generation of Schrodinger 
cats discussed in [25], but new interesting features stem 
from the spherical topology of the phase space and of 
the finite Hilbert space of the atomic sample. The coher- 
ent state superpositions are interesting not only per se as 
highly non-classical states, but were discussed as, e.g., a 
resource for building quantum logic gates [33] [34] . 

Let us assume that the Hamiltonian is -ff/ci of Eq. 
([7]) with f2i = fi| = flje and let the Hamiltonian be 



N-l 



for N even 
for N odd 



(23) 



During this time the basis states of J z and their superpo- 
sitions evolve as described in Appendix [B] For smoothly 
changing coefficients ha„ ~ ciAn+i of the superposition 
the total superposition contains pairs of states with one 
atom being either in the Rydberg state or in state \b). 
One can eliminate the Rydberg excitation by a pulse cou- 
pling the atomic states \b) and \ri) of duration tq with 



To = 



% for TV even 



n JC VN—i for iV odd 



(24) 



After that, the state contains almost exclusively contri- 
butions with no Rydberg excitation and with even An. 
The contribution of odd An states is suppressed to the 
order An 2 /N which is small for coherent states with An 
limited to An < \J~N j2. In particular, starting with a 
spin coherent state with N = 20, the resulting state has 
92% contribution of even An and 8% of odd An. Thus, 
one has created a state that is analogous to the even su- 
perposition of coherent states with oposite amplitudes. 

One can observe analogies of the most important fea- 
tures of the Schrodinger cat states familiar from optics: 
two separate peaks in one quadrature, interference fringes 
in the complementary quadrature, and "photon number" 
oscilations in the Fock representation. The results of our 
numerical simulation for N — 20 can be seen in Figs. 
[9-11 Although during the procedure the Rydberg state 
is populated with ~ 50% probability, after the last step 
the remaining population of the Rydberg state is only 
1.1%. In Fig. [9] one can see the oscillation in the atomic 
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FIG. 9: (Color online) Population of state \a) (wide blue 
bars) after the procedure of Schrodinger cat generation, with 
N = 20. The narrow red bars with tilted hatching show the 
population after rotation of the state by exp(i J x -w/2) show- 
ing the two separate peaks (live and dead cat). The narrow 
green bars with horizontal hatching show the population after 
rotation of the state by exp(i J v tt/2). 



population statistics (blue bars) analogous to the photon 
number oscillation in the even superposition of coherent 
states with opposite amplitudes. The red bars are the 
atomic population statistics after rotation of the result- 
ing state by exp(iJ x n/2). This is analogous to detecting 
an optical Schrodinger cat state by means of the homo- 
dyne scheme with the phase chosen to maximally sepa- 
rate the peaks. The green bars correspond to the state 
after rotation by exp(i J y it/2). In optics this would cor- 
respond to observing the state by homodyning such that 
the two peaks completely overlap and one detects the 
interference fringes. 

A Q-function of the resulting state can be seen in Fig. 
[TO] Only one component of the superposition is visible 
in the upper part of the figure, as the other one is on the 
opposite side of the sphere. One can see that the state 
is squeezed. To observe the squeezing by measuring the 
population of the atomic states, one can rotate the re- 
sulting state by exp(iJ y rf>) with a properly chosen angle 
4>. The resulting statistics is shown in Fig. [TTJ In this 
case one can see the green-bar statistics being narrower 
and the blue-bar statistics wider. In optical analogy this 
corresponds to a superposition of two phase-squeezed co- 
herent states with super Poissonian statistics. In this 
case the contribution of odd An of the resulting state 
dropped to 1.1%. 
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FIG. 10: (Color online) Q-function of the resulting state af- 
ter the Schrodinger cat generation procedure: two peaks on 
opposite sides of the sphere are produced. 
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FIG. 11: (Color online) Same as in Fig. Rafter the rotation 
of the state by exp(iJ !/ 0) where the angle 4> is optimized to 
have the statistics squeezed (narrow green bars with vertical 
hatching) or anti-squeezed (wide blue bars). 



VI. CONCLUSION 

We have proposed a scheme that takes advantage of 
the Rydberg blockade in a spin polarized atomic sam- 
ple to employ features of the JC model in a new envi- 
ronment. The finite atom number and the possibility 
to switch on and off the JC coupling between various 



atomic states makes the spectrum of interesting effects 
rather broad. In comparison to the standard JC model, 
one can squeeze coherent states symetrically, avoiding 
the "banana-like" deformation. By adiabatically chang- 
ing the Hamiltonian from the angular-momentum oper- 
ator to the double JC Hamiltonian, one can achieve full 
squeezing. With a proper choice of the JC Hamiltonian, 
one can prepare Schrodinger cat-like superpositions of 
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spin polarized states with opposite spins. These pro- the EU Integrated Project AQUTE. 
cedures can be useful, e.g., for precision time measure- 
ments, or for quantum information processing. 
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Appendix A: Parameters for the nonadiabatic transitions compensation 

Taking the norm of the vector (aiH ^ + a2H < f ( f ross ^ — Hi)\ip (t)} and expressing it as a function F(«i,a2) of (real) 
parameters 0*1,2 one gets 

F( ai ,a 2 ) = (M(<*iH% + a 2 H% cmss) - H 1 f\^) 

= (MH 2 M + aK^H^o) + al(i> \H% CTO ^ 2 \Tp Q ) - a 1 (^p \(H i f ( lH 1 + H.H^)^) 

-a% (-00 1 (-ff jc CTOSS) H i + ffi^c Cr ° SS) )|Vo) + ai a 2 (^\(H^ oss) H% + H%H^ oss) M ). (Al) 



Requiring for the extremum dF(cti,at2)/dcti — dF(ai,a2)/da>2 — one gets two equations for the unknowns 01,2 
with the solution 

0/ rr(y. cr ° ss ) 2 \/ tt(v) tt , tt tt{v)\ I rrfe, cross) „{y) , tt(v) tt(v, cross), , rriv, cross) „ „ rr(y, cross), 

a = 1 \ h .jc )\ h jc h i + h i h jc) " \ h jc h jc + h jc h jc )\ h jc h i + h i h jc ) ( A2 ) 

\ J C I \ J C 1 \ J C J C J C J C I 

a 2{H^ 2 ){H%^ Hl + Hl H%^) - (H^H% + H^H%^)(H^ Hl + H X H%) 

4<i4f ! ><ff& CTOS ^ 

where the mean value (...) is taken in state \tpo(t)}. 



Appendix B: Superposition generation for Schrodinger cat states 



Let us assume Hamiltonian Hjci of Eq. ([7| with fii = fi* = fljc being switched on for time r of Eq. (23). For 
even A, an initial state of the form | y + An, ^ — An, 0) evolves as (up to the second order in An) 



N A N A 
— + An, — - An, 



(-1) 



(JV+An)/2 



N . N 



An,0 



. An 2 ?r 
4A 



N A N 



An, 1 



for An even, and 



N A N A 
- + An, - - An. (I 



\ ^ (_ 1 )(JV+A«-l)/2 



in i 

4A 



N A N A 
- + An,y-An,0 



N A JV 

h An - 1, An, 1 

2 '2 



for An odd. For odd N, an initial state of the form 



N-l 

2 



An,%± 



An, 0) evolves as 



(Bl) 



(B2) 



(-1)^ 
for An even, and 



A- 1 
2 



An, 



A- 1 



An,0 ) +i 



4(A- 1) 



A- 1 « A — 1 

— 2 — + An, — An, 

A- 1 . A- 1 

h An - 1, An, 1 

2 '2 



(B3) 



Mr 



4(A-1) 



A- 1 



A N~l A 

An, An, 



A- 1 
2 

AT - 1 



A- 1 



An,0 



A- 1 
An - 1, An, 1 



(B4) 
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An 



for An odd. Thus, for even N a superposition of the form 

1*0) = J2 aA 

evolves into 

|* T> = (_l)"/2 / £ (_l)An/2 



N A N A n 
-+An,--An,0 



(B5) 



+ 



(An + 1) 2 7T 
47V 



OAn+1 



V An even 

N 



a An 



— + An, — - An, ) - ia An+1 



N A N A 

-+An,-- An -1,1 



- + An + 1, — - An - 1, ) - i± ^ a An+2 



N A iV A 

-+An + 1,--An-2,1 



and for odd N a superposition of the form 

An 

evolves into 



N-l A N-l , 

h An, An, 

2 2 



, (B6) 



(B7) 



£ (-D 



An/2 



v An even 



— JflAn+1 



iV-1 iV-1 \ (An+l) 2 7r 

— + An, — - An - 1, 1^) + 4(jV _ 1} °An + i x 



OAn 

iV - 1 



N-l a N-l A „ 

+ An, — An, 

N-l A 

-An + 1, — An -1,0 



.(An + 2) 2 7r 
'* 4(iV-l) 



OAn+2 x 



N-l A N — I 

— h An + 1, — An -2,1 



(B8) 
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